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ABSTRACT 


A  comparison  of  two  different  autopilot  configurations  and  their  effect  on  missile 
response  is  presented.  The  comparison  includes  the  steps  taken  in  determining  missile 
parameters  from  wind  tunnel  data  and  flight  condition  data.  The  missile  parameters  are 
coupled  with  two  different  autopilot  configurations  to  determine  any  significant  advan¬ 
tage  of  one  configuration  over  the  other.  Pole  placement  is  used  in  determining  required 
autopilot  feedback  and  feed  forward  gains.  Simulations  of  each  autopilot  are  conducted 
and  the  responses  are  compared. 
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I.  INTRODUCTION 


The  purpose  of  an  autopilot  in  a  missile  is  to  cause  the  missile  to  maintain  stability 
as  it  travels  along  its  flight  path.  Commonly  the  autopilot  provides  stability  about  the 
missile's  roll,  yaw  and  pitch  axes.  Autopilots  generally  use  displacement  gyros  with 
feedback  being  provided  by  electronic  sensors  which  detect  axial  displacement  due  to 
missile  rotation  about  one  or  more  of  the  above  mentioned  axes. 

Autopilots,  as  with  most  other  man-made  devices,  can  be  constructed  in  a  variety 
of  different  configurations.  Some  configurations  may  have  advantages  over  others  in 
terms  of  parameters  such  as  cost,  weight  or  effectiveness,  to  name  a  few. 

In  this  research  paper  two  pitch  axis  autopilot  configurations  are  compared.  The 
first  is  a  three  loop  autopilot  configuration,  employing  body  rate  feedback,  lateral  ac¬ 
celeration  feedback  and  an  additional  synthetic  stability  feedback  loop  as  depicted  in 
Figure  1. 


Figure  1.  Autopilot  Configuration  1 


The  second  is  a  two-loop  autopilot  configuration,  employing  proportional  plus  integral 
compensation  in  the  Iced  forward  path  with  body  rate  and  lateral  acceleration  feedback 
as  shown  in  figure  2. 


Figure  2.  Autopilot  Configuration  2 

In  these  two  figures  >/c  is  the  commanded  lateral  acceleration  in  the  pitch  plane,  tjM 
is  the  measured  lateral  acceleration  in  the  pitch  plane,  6  is  the  missile  control  surface 
deflection  angle  and  0  is  the  missile  pitch  angle  rate.  All  the  subscripted  K  variables  arc 
gains  which  arc  constant  at  each  design  point. 

The  purpose  of  the  comparison  of  the  two  autopilot  configurations  is  to  determine 
what  advantage,  if  any,  one  may  have  over  the  other  in  terms  of  missile  response. 
Comparisons  arc  made  of  the  two  autopilots  for  three  different  sets  of  flight  conditions. 
Each  set  of  flight  conditions  is  referred  to  as  a  design  point.  For  example,  design  point 
one  will  refer  to  the  first  set  of  flight  conditions.  These  flight  conditions,  or  design 
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points,  are  specified  in  Chapter  Two.  The  desired  autopilot  response  for  each  config¬ 
uration  is  given  as  approximately  10  radians  second  with  a  damping  coefficient  of  0.5. 

Pertinent  missile  parameters  are  derived  in  Chapter  Two.  using  given  wind  tunnel 
data  and  flight  condition  data.  Any  needed  parameters  which  cannot  be  derived  or 
otherwise  obtained  through  availible  information,  are  assumed  using  sound  engineering 
judgment. 

In  Chapter  Three  the  missile  and  autopilot  transfer  functions  are  developed.  The 
feedback  gains  necessary  to  meet  the  desired  autopilot  design  specifications  are  obtained 
allowing  development  of  simulation  studies. 

Simulations  are  developed  and  conducted  using  Dynamic  Simulation  Language 
(DSL)  in  Chapter  Four  and  the  results  are  analyzed. 


3 


II.  DEVELOPMENT  OF  MISSILE  PARAMETERS 


A.  INTRODUCTION 

In  order  to  compare  the  effects  that  the  two  autopilot  configurations  have  on  missile 
performance,  determination  of  certain  missile  parameters  must  first  be  made.  Specif¬ 
ically.  the  transfer  functions  0{s)lci{s)  and  (s)/S(s),  which  are  the  missile  dynamics 
blocks  in  the  autopilot  block  diagrams  of  Figure  1  on  page  1  and  Figure  2  on  page  2, 
must  be  obtained. 

In  order  to  determine  these  transfer  functions,  wind  tunnel  data  and  three  sets  of 
flight  condition  data  are  provided  later  in  this  chapter.  The  procedures  for  processing 
this  data  are  also  presented  later  in  this  chapter. 

Prior  to  developing  the  missile  parameters,  an  explanation  of  the  notation,  termi¬ 
nology,  and  coordinate  system  used  is  in  order.  [Ref.  1:  p.  132]  Refer  to  Figure  3. 
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The  notation  in  Figure  3  represents  the  following: 

CG  is  the  center  of  gravity  of  the  missile. 

CP  is  the  aerodynamic  center  of  pressure  of  the  missile. 
l'M  is  the  velocity  of  the  missile. 

SRL  is  the  spatial  reference  line,  a  reference  fixed  in  space. 

Fl  is  the  lift  force.  It  is  perpendicular  to  the  velocity  vector. 

W  is  the  weight  of  the  missile. 

MP  is  the  pitching  moment.  It  is  taken  as  positive  in  the  counter-clockwise  direction. 

a  is  the  angle  of  attack.  It  is  the  angle  taken  from  the  velocity  vector  to  the  longi¬ 
tudinal  missile  axis. 

6  is  the  angle  between  the  SRL  and  the  longitudinal  missile  axis. 

d  is  the  fin  deflection  angle.  It  is  taken  from  the  longitudinal  missile  axis  to  the  control 
fin  axis. 

P  is  the  angle  between  the  gravity  vector  and  the  SRL. 

Other  notation  and  terminology  will  be  discussed  as  it  is  introduced.  It  is  important  to 
remember  that  this  problem  considers  only  motion  in  the  pitch  plane  but  that  similar 
methods  may  also  apply  to  yaw  motion. 

B.  FLIGHT  CONDITION  DATA 

The  three  sets  of  flight  condition  data  given  in  Table  1  on  page  6  establish  the  de¬ 
sign  points  around  which  the  autopiloted  missile  will  be  examined.  The  flight  condition 
data  will  be  used,  in  conjunction  with  wind  tunnel  data  presented  in  the  next  section,  to 
establish  the  missile  parameters  and  autopilot  gains  necessary  to  develop  and  conduct 
simulations  of  the  autopiloted  missile. 

Some  manipulation  of  the  flight  condition  data  is  required  in  order  for  it  to  be  used. 
These  conversions  and  calculations  will  be  made  as  necessary  and  explained  when  used. 
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Table  1.  FLIGHT  CONDITIONS  AT  THREE  DE¬ 
SIGN  POINTS 


DESIGN  POINT  1 


T  =  Time  (referenced  to  launch)  =  9.5  s 
Alt  =  Altitude  =  4.993  km 
V  =  Velocity  =  1112.4  ms 
M  =  Mass  =  450  kg 
q  =  Dynamic  Pressure  =  455642  N  mJ 
I  =  Moment  oflnertia  =  720  kg  m2 
Power  On 

CP  is  265  cm  aft  of  nose 


DESIGN  POINT  2 


T  =  18.5  s 
Alt  =  13.54  km 
V  =  1086.4  ms 
M  =  370  ka 
q  =  146S29  N  m2 
I  =  693  kg  mJ 
Power  Oil 

CP  is  265  cm  aft  of  nose 


DESIGN  POINT  3 


T  =  22.5  s 
Alt  =  17.82  km 
V  =  934.5  ms 
M  =  365  kg 
q  =  54637  N  m3 
I  =  6S7  kg  m: 

Power  Off 

CP  is  266  cm  aft  of  nose 


C.  WIND  TUNNEL  DATA 

The  wind  tunnel  data  presented  in  this  section  will  be  instrumental  in  determining 
missile  parameters.  The  data,  shown  in  Table  10  on  page  37  through  Table  13  on  page 
40,  are  extracted  from  the  complete  wind  tunnel  data  set  for  the  missile  used  in  this  re¬ 
search  paper.  These  data,  as  well  as  the  three  sets  of  flight  condition  data  contained  in 
Table  1,  were  provided  by  the  U.S.  Army's  Missile  and  Space  Intelligence  Center 
(MSIC). 

The  first  two  sets  of  wind  tunnel  data  are  the  axial  force  coefficients.  These  are 
given  in  Table  10  and  Table  1 1  for  the  Power  On  and  Power  Off  cases. 
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The  third  set  of  wind  tunnel  data  is  the  normal  force  coefficients.  These  data  in 
Table  12  on  page  39  will  be  used,  along  with  the  axial  force  coefficient  data  to  calculate 
the  lift  and  drag  of  the  missile  at  the  design  points. 

The  final  set  of  wind  tunnel  data  is  the  pitching  moment  coefficients.  These  data, 
in  Table  13  on  page  40  will  be  used  to  determine  the  pitching  dynamics  of  the  missile. 

D.  WIND  TUNNEL  DATA  INTERPOLATION 

To  enter  the  wind  tunnel  data  at  the  desired  design  point,  the  Mach  number  must 
be  calculated  from  the  flight  condition  data.  Using  given  missile  velocity  and  altitude, 
with  temperature  obtained  from  gas  tables  and  graphs  found  in  [Refs.  2,3],  the  Mach 
number  for  each  design  point  is  calculated  by  gas  law  equation  techniques.  These  cal¬ 
culations  are  not  included  here.  The  Mach  number  at  each  design  point  is  given  in 
Table  2. 


Table  2.  MACH  NUMBERS  AT  DESIGN  POINTS 


DESIGN  POINT 

MACH  NUMBER 

1 

3.4S 

2 

3.68 

> 

3.17 

The  availability  of  the  Mach  number  at  each  design  point  allows  entry  of  the  wind 
tunnel  data  and  interpolation  of  the  various  coefficients  for  each  design  point.  The  in¬ 
terpolation  of  the  wind  tunnel  data  at  each  design  point  is  accomplished  using  a  simple 
linear  interpolation  scheme.  The  interpolation  provides  the  coefficients  at  each  design 
point.  These  interpolated  coefficients  are  given  in  Table  14  on  page  41  through 
Table  17  on  page  43. 

E.  CALCULATION  OF  LIFT  AND  DRAG  COEFFICIENTS 

Once  the  axial  and  normal  force  coefficients  at  the  design  points  have  been  deter¬ 
mined  the  lift  and  drag  coefficients  can  be  calculated.  The  lift  force,  FL,  and  drag  force, 
F0  ,  are  related  to  the  normal  force,  FN,  and  the  axial  force,  FA,  by  the  trigonometric 
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relationships 


FL  =  FyCosa  —  FAs  inoc 


(2.1) 


and 


FL)  =  /\vsina  +  F/(cosa.  (2.2) 

These  relationships  are  depicted  in  Figure  4  below. 


Since  the  coefficients  are  dimensionless  representations  of  the  forces,  the  same  re 
lationships  apply.  Specifically, 


and 


C,  =  CAcosa  -  Qsina  (2.3) 


CD  =  CAsina  +  CA  cosa  (2.4) 
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where  CL  is  the  lift  coefficient,  CD  is  the  drag  coefficient,  CA  is  the  axial  force  coefficient, 
and  C.  is  the  normal  force  coefficient.  [Ref.  1:  p.  134] 

The  application  of  these  equations  to  the  interpolated  wind  tunnel  data  in  Appendix 
B  renders  lift  and  drag  coefficients  at  the  design  points.  These  lift  and  drag  coefficients 
are  tabulated  in  Appendix  C. 

F.  ANALYSIS  OF  WIND  TUNNEL  DATA 

There  are  four  missile  parameters  which  are  of  interest  in  the  development  of  the 
missile  dynamics  transfer  functions.  These  parameters  are:  Fu,  which  is  the  lift  force  per 
angle  of  attack,  FLi.  which  is  the  lift  force  per  fin  deflection  angle,  MP%,  which  is  the 
moment  about  the  center  of  gravity  per  angle  of  attack,  and  Mp(,  which  is  the  moment 
about  the  center  of  gravity  per  fin  deflection  angle. 

These  parameters  can  be  determined  from  their  corresponding  coefficients  and  a  few 
physical  dimensions  of  the  missile.  The  corresponding  coefficients  are  discussed  in  the 
paragraphs  which  follow. 

CL,  is  the  lift  coefficient  per  angle  of  attack.  It  is  determined  by  graphing  the  lift 
coefficients  against  the  angle  of  attack  and  approximating  the  slope  of  the  nearly  linear 
graph. 

Cu  is  the  lift  coefficient  per  fin  deflection  angle.  It  is  determined  by  graphing  the  lift 
coefficients  against  the  fin  deflection  angle  and  approximating  the  slope. 

Similarly,  CM%  and  CM,  are  the  pitching  moment  coefficients  per  angle  of  attack  and 
fin  deflection  angle,  respectively.  They  too  are  obtained  by  graphing  the  coefficients 
against  the  angles  and  approximating  the  slope  of  the  resulting  graph. 

The  graphs  of  all  the  coefficients  versus  corresponding  angles  are  contained  in  Ap¬ 
pendix  D.  The  graphs  are  made  for  each  set  of  coefficients  at  each  of  the  three  design 
points.  The  coefficients  Cu,  CLi,  CM% ,  and  CMi  are  tabulated  in  Table  3  on  page  10. 
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Table  3.  COEFFICIENTS  DUE  TO  ANGLE  OF  ATTACK  AND 
FIN  DEFLECTION  ANGLE 


DESIGN 
POINT  1 

DESIGN 
POINT  2 

DESIGN 
POINT  3 

CLx 

0.27 

0.25 

0.28 

CL6 

0.03 

0.02 

0.03 

Q/a 

-0.15 

-0.12 

-0.20 

C\I6 

-0.17 

-0.15 

-0.22 

G.  CALCULATION  OF  MISSILE  PARAMETERS 

As  mentioned  in  the  preceding  section,  the  parameters  of  interest  are  Fu,  Fu,  MPt, 
and  3/r,s.  With  knowledge  of  the  lift  and  moment  coefficients  previously  obtained,  and 
knowledge  of  some  basic  physical  dimensions  of  the  missile,  these  parameters  are  deter¬ 
mined.  The  necessary  missile  dimensions  are  shown  in  Figure  5.  These  missile  dimen¬ 
sions  were  provided  by  the  Missile  and  Space  Intelligence  Center. 


Figure  S.  Missile  Dimensions 


ft 
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The  equations  used  to  determine  the  missile  parameters  from  the  coefficients  are: 

r^C^qA,  (2.5) 

F L6  =  C-uflA,  (2.6) 

Up*=CujAc,  (2.7) 

and 

MP6  =  CM6qAc  (2.8) 

where  q  is  the  dynamic  pressure,  provided  in  the  flight  conditions,  A  is  the  missile  char¬ 
acteristic  area,  and  c  is  the  characteristic  length.  The  characteristic  area,  A,  is  taken  as 

A  =  n-£-  (2.9) 

where  d  is  the  missile  diameter.  38  cm,  as  given  in  Figure  19  above.  The  characteristic 
length  is  taken  as  the  missile  diameter.  (Ref.  1:  p.  134] 

Applying  the  above  equations  to  the  known  coefficients  produces  the  missile  pa¬ 
rameters  at  the  three  design  points.  Use  of  these  parameters  in  later  chapters  requires 
that  they  be  given  in  force  or  moment  per  radian  rather  than  per  degree.  This  conver¬ 
sion  is  made  mulitplying  each  parameter  by  ISO /n  .  These  converted  parameters  are 
tabulated  in  Table  4  on  page  12. 

Knowledge  of  these  missile  parameters  will  enable  the  derivation  of  the  missile  dy¬ 
namics  transfer  functions  in  the  following  chapter. 
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Table  4.  MISSILE  PARAMETERS  AT  DESIGN  POINTS 


PARAMETER 

DESIGN 
POINT  1 

DESIGN  ; 
POINT  2 

DESIGN 
POINT  3 

799276  N 

238350  N 

99695  N 

Ft* 

SSSOS  N 

19080  N 

10657  N 

MPy 

-16S736  Nm 

-43545  Nm 

-269S6  Nm 

Ups 

-19136S  Nm 

-54431  Nm 

-29679  Nm 

III.  DETERMINATION  OF  TRANSFER  FUNCTIONS 


A.  OVERVIEW 

In  this  chapter,  the  missile  parameters  obtained  in  the  previous  chapter  are  applied 
to  the  equations  of  motion  which  govern  the  flight  of  the  missile.  From  this,  the  missile 
dynamics  transfer  functions.  8(s)IS(s)  and  t]M(s)IS(s),  are  derived.  By  then  inserting  these 
two  missile  dynamics  transfer  functions  into  the  two  autopilot  transfer  functions,  the 
complete  autopilot  transfer  functions  can  be  expressed  in  terms  of  missile  parameters 
and  feedback  gains. 

B.  MISSILE  DYNAMICS  TRANSFER  FUNCTIONS 

Consider  the  forces  and  moments  acting  on  a  missile  while  in  flight.  These  are  de¬ 
picted  in  Figure  3  on  page  4.  Summing  the  forces  in  the  transverse  direction,  that  is 
perpendicular  to  the  velocity  vector,  yields 


y \ft  =  t]um  —  Fl~  IV  sin(/?  +  6  -  a). 


(3.1) 


Since 


fl  =  flFj  +  FiA&  +  «)» 


(3.2) 


the  substitution  is  made  yielding 

*}\tm  =  FU°-  +  FiA*>  +  «)  -  ,r  sin(/J  +  8  -  a).  (3.3) 

Dividing  through  by  the  mass,  m,  yields  the  lateral  acceleration,  . 

1*/“'Tnmec+-jn~(t  +  *) — sin  (fi  +  0-«)  (3-4) 

The  weight  of  the  missile  is  neglected  at  this  point  for  simplicity.  The  missile's 
weight  would  affect  the  lateral  acceleration  of  the  missile  most  when  the  angle 
(/?  +  8  —  a)  is  90  or  270  degrees.  At  this  angle  the  contribution  to  lateral  acceleration 
would  be  plus  or  minus  1  G.  This  effect  would  be  offset  by  fin  deflection  and  angle  of 
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attack  trim  angles  of  much  less  than  one  degree  at  any  of  the  three  design  points. 
Neglecting  missile  weight  leaves  the  approximation 


V\i— 


1  Lr 

m 


I  ^LS  ,  -  v  +  FLfi)  (FL6) 

a  +  —777—  (d  +  a)  =  a - - - +  d  ■ 


m 


m 


(3.5) 


which  is  used  to  determine  rjM  at  each  design  point  in  terms  of  a.  and  <5. 

In  a  similar  manner,  summing  moments  about  the  center  of  gravity  of  the  missile 
with  the  counter-clockwise  direction  being  referenced  as  the  positive  direction  yields 


Va 1cg  =  16  =  MP  (3.6) 

where  I  is  the  missile  s  moment  of  inertia,  9  is  the  missile's  angular  acceleration  in  the 
pitch  plane,  and  MP  is  the  pitching  moment.  The  pitching  moment,  Mp,  is  given  by  the 
equation 


Mp  —  aMp,  +  9Mpe  +  aA/p,  +  (a  +  d)MP6.  (3.7) 

Since  a  0  and  a  ss  6  the  substitution  is  made  leaving 

MP  =  a{MPx  +  Mpe)  +  a(.Vp2  -I-  MPS)  +  S(MP6).  (3.S) 


Combining  and  rearranging  the  above  equations  yields 

(•^4>  +  -Wpe)  .  (Mp*  +  Mps)  Mps 
* - ] - a - j - a  =  — j —  5.  (3.9) 

In  Equation  (3.9),  the  term 

Wh  + 

I 

is  the  aerodynamic  damping  term.  It  is  beyond  the  scope  of  this  research  paper  to 
properly  obtain  this  term.  It  is,  therefore,  assumed  that  the  missile  is  designed  with  a 
damping  coefficient  of  (  =  0.5.  The  aerodynamic  damping  term  will  be  different  nu¬ 
merically  at  each  of  the  three  design  points  and  will  be  referred  to  hereafter  as  D.  It  will 
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be  calculated  later  in  this  chapter.  [Ref.  4] 
Substituting  D  into  Equation  (3.9)  yields 


a  —  Dy.  — 


(Ml'y  +  '//>t0 


.1/ 


I 


■  a  =  ■ 


(3.10) 


Using  LaPlacc  transforms  and  assuming  zero  initial  conditions  this  equation  transforms 
into 


Rearranging,  this  becomes 


Ups 

oefs) _ / _ 

(5(s)  ^  (.!//»,  +  Mp.O 

s'  —  Ds - : - — 


Multiplying  both  sides  by  s  renders 


sa(.s)  _ 5  / _ «(•?) 

d(s)  2  \Mp,  +  Mpt) 

s  -  Ds - j - 

Since  a  =  6, 

Mrs 

a(s)  6(s) _ 5  I _ 

d(s)  d(s)  2  (MPx  +  Mp6) 

s  —  Ds - * — ; - 


(3.11) 


(3.12) 


(3.13) 


(3.14) 


This  transfer  function  is  depicted  in  the  signal  flow  graph  of  Figure  6  on  page  16. 
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Figure  6.  Pitch  Transfer  Function  Signal  Flow  Graph 

The  aerodynamic  damping  term,  D,  is  calculated  for  each  of  the  three  design  points. 
The  denominator  of  the  0(s)IS{s)  transfer  function  is  the  characteristic  equation  of  that 
system.  It  is  in  the  quadratic  form  of  s 1  •+  s2twn  +  Equating  coefficients  yields 

-D  =  2(wn  (3.15) 


and 

I  Wp*  +  n  m 

- j j - •  (3-16) 

Using  the  missile  parameters  obtained  in  the  previous  chapter  and  the  assumed  value 
of  0.5  for  £,  the  damping  term,  D,  is  calculated.  The  values  for  D  are  given  in  Table  5 
on  page  17.  (Ref.  5:  pp.  106-107] 
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Table  5.  AERODYNAMIC  DAMPING  COEFFI¬ 
CIENTS 


DESIGN  POINT 

DAMPING  COEFFICIENT 

1 

-22.4 

*) 

-1 1.9 

3 

-9.1 

The  other  missile  dynamics  transfer  function  needed  is  the  lateral  acceleration 
transfer  function.  VM(s)!S{s).  Combining  the  signal  How  graph  of  Figure  6  with 
Equation  (3.5)  results  in  the  signal  How  graph  of  Figure  7. 


e 


Figure  7.  Acceleration  Transfer  Function  Signal  Flow  Graph 


The  transfer  function  associated  with  this  signal  flow  graph  is  given  as 

f)  _  I  "i _  rLfi 

2  „  (-Vh  +  AW  ' 

s  —  sD - ; - 


(3.17) 
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The  transfer  functions  for  the  missile  dynamics  are  thus  obtained.  They  will  be  used 
in  the  following  section  to  determine  the  autopilot  transfer  functions. 

C.  AUTOPILOT  TRANSFER  FUNCTIONS 

By  inserting  the  two  missile  dynamics  transfer  functions  into  each  of  the  two 
autopilot  configuration  block  diagrams  presented  in  Chapter  One,  and  performing  a  few 
basic  block  diagram  manipulations,  the  closed  loop  autopilot  transfer  functions  are  ob¬ 
tained.  In  both  cases,  the  desired  autopilot  transfer  function  is  >iM(s)l>]c(s).  For  Autopilot 
Configuration  One  the  transfer  function  is 


V\f{s) 

'/cO) 


d(s) 


7  .  .  Q(s)  6{s) 

S  +  s(kF  +  kBR  *7—7  )  +  KST&8R  1 


d(s) 


+  KA 


(3.18) 


Expanding  the  missile  dynamics  transfer  functions  in  the  autopilot  transfer  function  and 
performing  extensive  algebraic  manipulations  yields 


G[(s) 


V  \f(s) 
>/cW 


n0]s  +nns  +  n 


'2! 


5  +  ]  S  +  t/j  j  S  +  (i^  j  S  +  t/3  j 


(3-19) 


where 


«oi  - 


=  -D 


1  L6 

m 


Mp*  IFl*  +  Fu>) 

j  ni 


F 16  (Mp*  +  Mpi) 

in  1 


^01  =  kf—  d  , 


M 


4,  =  KBR-p-+KA-^-KFD- 


LS 


(Mp3  +  Mp6) 


—  KSTKBR  ' 


Mpi  „  (.V/Pa  +  MP6)  v  nt-u 


and 


IS 


,  „  ,  3//*  (^Za  +  ^uO  S  (-'/fj  +  -W/\0  , 

aii~hA\  j  ,n  ~  ,n  j  )■ 

The  transfer  function  for  Autopilot  Configuration  Two  is  obtained  in  a  similar 
manner.  For  the  second  configuration 


>lds) 


{K2s  +  A',) 


>7uU) 

h(s) 


2  ,  ,  rr 

5  +  s{kr+k™~) 


)+KA(K2s+k\) 


Vm(s) 

<3(5j 


This  expands  into 


G2(s)  = 


>lds) 


h$2s  4-  W]25*  +  "22s  +  " 


32 


*02J 


d\2s  +  d22s  +  c/22 


(3.20) 


(3.21) 


where 


Wq2  —  ^2 


r 


_Z<5_ 

m 


i'  kl-6  i'  nJlLL 

"\2  “i  in  h2u  /u 


,,  ,  Sfp6  (Fu  +  Fls)  Fls  (-'/ft,  +  3;w)  % 

"22  =  hi— - m - m - / - )~hD 


in 


_  t'  ( ‘On  L Li + [Mfi + Mpd  s 

w32  ^  m  I  )  ’ 


I 


ta-Kf-D+K^-!#- 


Mp6  Mp  +  Mpfi  Fix  Fix 

dn  =  K,x-T- - ,  --KFD-KJK2D-g-+KilKI-g-  . 


/ 


di2  =  KAK2{ 


x ,  Up6  {Fu  +  f^)  Fu  (Mfa  +  Mn) 


m 


m 


)-KAKlD-^--KF 


Fls  „  (MPa  +  Mw) 


and 


19 


rfjj  -  *>',( 


-»«  +  fm)  fl,'  W,  + 

/  in  m  i  >  ‘ 

These  closed  loop  autopilot  transfer  functions  will  be  used  in  the  next  section  for 
system  pole  placement  and  solution  of  the  feedback  and  feed  forward  gains. 

D.  POLE  PLACEMENT  AND  GAIN  DETERMINATION 

As  stated  in  Chapter  One  of  this  research  paper,  the  desired  autopilot  response  for 
both  configurations  is  con  =  10  with  a  damping  coefficient  £  =  0.5  .  These  design  spec¬ 
ifications  dictate  a  complex  pair  of  poles  in  the  characteristic  equation  whose  real  part 
is  -5.0.  Since  the  autopilots'  characteristic  equations  are  fourth  order,  two  more  poles 
are  required  to  put  them  in  the  form 

CE  =  (s 2  +  s2£w„  +  cj^)(s  +  /ii)(s  +  />2) .  (3.22) 

In  order  for  the  complex  pair  placed  by  design  specifications  to  be  dominant,  the  other 
two  poles  must  lie  to  the  left  of  them  in  the  s-plane.  The  two  additional  poles  aie  se¬ 
lected  for  simplicity  as  real  and  equal  and  placed  at  -20.0.  The  characteristic  equation 
for  both  autopilot  configurations  is 

CE=(s  +  20  )(s  +  20)(s2  +  10s  +  100) 

=  s4  +  50s3  +  900s2  +  SOOOs  +  40000  .  (3.23) 

The  denominators  of  the  closed  loop  transfer  functions,  which  are  given  by 
Equations  (3.19)  and  (3.21).  are  the  autopilot  characteristic  equations.  The  transfer 
function  gains  are  solved  for  by  first  equating  the  coefficients  of  characteristic  equation 
of  above  to  those  of  Equations  (3.19)  and  (3.21),  then  solving  simultaneous  equations 
for  the  gains.  The  calculated  gains  for  Autopilot  Configuration  One  and  Autopilot 
Configuration  Two  are  tabulated  in  Table  6  and  Table  7,  respectively.  The  determi¬ 
nation  of  the  autopilot  gains  makes  it  possible  to  simulate  the  response  of  the  missile. 
These  simulations  will  be  conducted  in  the  next  chapter.  (Refs.  6,7] 
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Table  6.  AUTOPILOT  GAINS  FOR  CONFIGURATION  ONE 


GAIN 

DESIGN 
POINT  1 

DESIGN 
POINT  2 

DESIGN 
POINT  3 

-0.094 

-0.S45 

NO 

o 

1 

Kf 

27.6 

3S.1 

40.9 

k-ST 

26.96 

8.9S 

10.18 

^ BR 

0.752 

-4.44 

-12.37 

Table  7.  AUTOPILOT  GAINS  FOR  CONFIGURATION  TWO 


GAIN 

DESIGN 
POINT  1 

DESIGN 
POINT  2 

DESIGN 
POINT  3 

A^A, 

-0.094 

-0.845 

-3.06 

A,  A, 

0.004 

-0.057 

-0.352 

20.1 

41.05 

51.17 

k-BR 

-0.513 

-7.62 

-12.S3 

It  is  interesting  to  note  that  the  values  of  A',  and  K2  for  Configuration  Two,  depend 
on  the  selected  value  of  A*.  For  purposes  of  simulation  in  the  next  chapter,  KA  is  chosen 
as  unity.  This  results  in  A* A,  being  equal  to  A':  and  KAK2  being  equal  to  K2  for  the  sim¬ 
ulation. 
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IV.  AUTOPILOT  SIMULATION  AND  ANALYSIS 


A.  AUTOPILOT  SIMULATION 

In  order  to  compare  the  responses  of  the  two  autopilot  configurations,  a  simulation 
program  was  developed  using  IBM's  Dynamic  Simulation  Language.  [Ref.  8] 

The  program  simulates  the  step  response  of  the  autopilot  transfer  functions  for  both 
configurations,  at  all  three  design  points.  The  simulation  was  developed  from  Equations 
(3.19)  and  (3.21),  which  describe  the  two  autopilot  transfer  functions,  and  uses  the 
feedback  and  feed  forward  gains  tabulated  in  Table  6  and  Table  7  on  page  21.  The 
program  also  incorporates  a  saturating  limiter  which  limits  the  deflection  of  the  control 
fins  to  ±  20°.  A  copy  of  the  program  listing  is  at  Appendix  E. 

The  program  was  run  on  the  IBM  3360  mainframe  computer  system  at  the  Naval 
Postgraduate  School  utilizing  a  Tektronix  618  Monitor  and  Tektronix  4631  Hard  Copy 
Unit  for  graphic  output. 

This  chapter  contains  the  graphical  results  of  these  simulations.  Figure  8  on  page 
23  through  Figure  10  on  page  25  depict  the  step  response  of  Autopilot  Configuration 
One  at  each  of  the  three  design  points.  Figure  1 1  on  page  26  shows  the  parametric 
display  of  the  step  response  of  Autopilot  Configuration  One  at  all  three  design  points. 
Figure  12  on  page  27  through  Figure  14  on  page  29  depict  the  step  response  of 
Autopilot  Configuration  Two  at  each  of  the  three  design  points.  Figure  15  on  page  30 
displays  parametrically  the  step  response  of  Autopilot  Configuration  Two  at  all  three 
design  points.  The  parametric  comparisons  of  the  two  autopilot  configurations  at  each 
design  point  are  shown  in  Figure  16  on  page  31  through  Figure  18  on  page  33.  The 
next  section  contains  discussion  and  analysis  of  the  graphical  output. 
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■  • 


Figure  14.  Step  Response  of  Autopilot  Configuration  Two  at  Design  Point  Three 


Figure  15 


Parametric  Comparison  of  the  Step  Response  of  Autopilot  Configuration 
Two  at  Three  Design  Points 
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STEP  RESPONSE.  CONFIGURATIONS  1  &  2.  DESIGN  POINT  2 


of  Autopilot  Configure- 


B.  ANALYSIS  OF  SIMULATION  RESULTS 


In  the  analysis  of  the  simulation  results,  there  are  two  characteristics  of  the  step  re¬ 
sponse  which  are  of  interest.  They  are  the  rise  time.  ir ,  and  the  settling  time,  i  . 

The  rise  time  is  simply  defined  as  the  time  required  for  the  systems  output  response 
to  go  from  10  percent  to  90  percent  of  its  final  value.  This  can  be  measured  directly  oft 
of  the  graphical  simulation  results.  [Ref.  9:  p.  40] 

The  settling  time,  r: ,  is  defined  as  four  time  constants,  or 


(4.1) 


The  value  of  the  maximum  overshoot.  Mpl ,  and  the  time  at  which  the  maximum  over¬ 
shoot  occurs,  ip  ,  can  be  measured  directly  off  the  graphical  output.  The  maximum 
overshoot  is  defined  as 


Mpt  =  1  +  exp( 


(4.2) 


and  the  time  at  which  it  occurs 


"  T  -2 

Un\  1  ~  * 

Measurements  of  Mpt  and  ip  are  made  from  the  graphical  output  of  Figure  8  on  page 
23  through  Figure  10  on  page  25  for  Configuration  One  and  from  Figure  12  on  page 
27through  Figure  14  on  page  29  for  Configuration  Two.  Once  Mp,  and  ip  are  measured, 
£  and  o>„  are  obtained  by  solving  simultaneously  Equations  (4.2)  and  (4.3).  The  values 
for  £  and  tu„  are  then  used  in  Equation  (4.1)  to  obtain  the  settling  time,  The  response 
characteristics  Mpt,  tp,  £,  co„ ,  t,  and  /,  are  tabulated  in  Table  8  and  Table  9  on  page  35. 
[Ref.  5:  pp.  108-109] 
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Table  8.  RESPONSE  CHARACTERISTICS  FOR  AUTOPILOT 
CONFIGURATION  ONE 


DESIGN 
POINT  1 

DESIGN 
POINT  2 

DESIGN 
POINT  3 

1.12 

1.13 

1.07 

lp 

0.51  sec 

0.51  sec 

0.57  sec 

0.56 

0.54 

0.65 

(°n 

7.47  rad'sec 

7.36  rad  sec 

7.24  rad'sec 

‘s 

0.96  sec 

1.01  sec 

0.S5  sec 

' r 

0.21  sec 

0.19  sec 

0.23  sec 

Table  9.  RESPONSE  CHARACTERISTICS  FOR  AUTOPILOT 
CONFIGURATION  TWO 


DESIGN 
POINT  1 

DESIGN 
POINT  2 

DESIGN 
POINT  3 

Up, 

1.13 

1.19 

1.24 

h 

0.63  sec 

0.52  sec 

0.43  sec 

Y 

0.54 

0.47 

0.41 

<*>n 

5.95  rad  sec 

6.84  rad  sec 

8.03  rad,  sec 

h 

1.24  sec 

1.24  sec 

1.20  sec 

tr 

0.23  sec 

0.19  sec 

0.14  sec 
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V.  CONCLUSIONS 


The  step  responses  of  the  two  autopilot  configurations  were  compared  at  three  dif¬ 
ferent  design  points.  Both  configurations  met  the  desired  response  specifications  of 
(2:0.5  and  a>„2:10.  The  actual  values  of  (  and  &>„  are  given  in  Table  8  and  Table  9  on 
page  35  for  Configurations  One  and  Two  respectively. 

At  Design  Point  One,  Configuration  One's  observed  rise  time  is  nine  percent  less 
than  that  of  Configuration  Two,  and  Configuration  One's  observed  settling  time  is 
twenty-three  percent  less  than  that  of  Configuration  Two.  The  two  configurations  are 
compared  graphically  in  Figure  16  on  page  31. 

At  Design  Point  Two.  Configuration  One's  observed  rise  time  is  equal  to  that  of 
Configuration  Two,  and  Configuration  One's  observed  settling  time  is  nineteen  percent 
less  than  that  of  Configuration  Two.  The  graphical  comparison  is  shown  in  Figure  17 
on  page  32. 

At  Design  Point  Three,  Configuration  One's  observed  rise  time  is  sixty-four  percent 
more  than  that  of  Configuration  Two,  and  Configuration  One's  observed  settling  time 
is  thirty  percent  less  than  that  of  Configuration  Two.  This  comparison  is  shown  graph¬ 
ically  in  Figure  18  on  page  33. 

The  apparent  trend  is  that  Configuration  One’s  rise  time  increases  slightly  as  the 
flight  of  the  missile  continues  but  that  its  settling  time  decreases  slightly.  This  trend  is 
graphically  depicted  in  Figure  11  on  page  26.  Configuration  Two  s  rise  time  decreases 
significantly  as  the  flight  of  the  missile  continues  but  its  settling  time  remains  nearly 
constant.  This  trend  is  shown  graphically  in  Figure  15  on  page  30. 

The  shorter  rise  time,  later  in  flight,  of  Configuration  Two,  is  viewed  as  a  significant 
performance  advantage  over  Configuration  One.  It  is  clear  that  a  shorter  rise  time 
means  a  faster  response.  This  faster  response  is  desirable  late  in  flight  because  it  is 
during  this  later  part  of  the  missile's  flight  that  the  missile  tvould  be  in  its  final  phase  of 
closing  on  and  tracking  a  maneuvering  target.  For  this  reason,  faster  response,  later  in 
flight,  is  a  significant  advantage. 
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APPENDIX  A.  ORIGINAL  WIND  TUNNEL  DATA 


Table  10.  AXIAL  FORCE  COEFFICIENTS  (POWER  ON) 


ALT  (KM) 
(kmi 

MACH  2.60 

MACH  3.20 

MACH  3. SO 

0 

0.324 

0.25S 

3.0 

0.326 

MiM 

0.260 

0 

5.0 

- 

0.292 

0.263 

10.0 

0.305 

0.278 

20.0 

0.35S 

0.336 

0.322 

30.0 

0.24S 

0.323 

0.422 

0 

HI 

0.295 

WKBSM 

3.0 

0.296 

Mis  II  • 

8 

5.0 

mmm 

0.299 

10.0 

0.348 

0.312 

-w%m 

20.0 

0.367 

0.343 

WmmM 

30.0 

0.257 

0.330 

0.427 

0 

MEBM 

3.0 

■EL'. 

■ 

WMEm 

16 

5.0 

0.357 

0.316 

10.0 

0.367 

0.329 

0.297 

20.0 

0.3S7 

0.360 

0.342 

30.0 

0.276 

0.347 

0.441 

0 

mmm 

mmm 

■Sim 

3.0 

24 

5.0 

msEm 

0.340 

0.304 

10.0 

0.395 

0.353 

0.319 

20.0 

0.415 

0.3S4 

0.363 

30.0 

0.304 

0.371 

0.462 
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Table  12.  NORMAL  FORCE  COEFFICIENTS 


DELTA 

ALPHA 

( decrees) 

MACH  2.60 

0 

■BUSH 

-0.999 

-0.525 

3.0 

-0.392 

-0.015 

-30 

5.0 

0.013 

0.353 

10.0 

1.363 

1.594 

20.0 

5.004 

4.877 

4.739 

30.0 

9.420 

8.992 

8.391 

0 

-0.802 

-0.660 

-0.423 

3.0 

-0.138 

-0.041 

0.102 

-20 

5.0 

0.291 

0.36S 

0.479 

10.0 

1.699 

1.707 

1.743 

20.0 

5.358 

5.163 

4.926 

30.0 

9.758 

9.259 

8.604 

0 

-0.385 

-0.305 

-0.236 

3.0 

0.272 

0.294 

0.300 

-10 

5.0 

0.6S9 

0.692 

0.683 

10.0 

2.076 

2.008 

1.959 

20.0 

5.723 

5.451 

5.159 

30.0 

10.156 

9.577 

8.839 

0 

0.000 

o.ooo 

0.000 

3.0 

0.652 

0.594 

0.539 

0 

5.0 

1.067 

0.9S9 

0.923 

10.0 

2.457 

2.309 

2.201 

20.0 

6.142 

5.7SS 

5.391 

30.0 

10.625 

10.004 

9.043 

0 

0.3S5 

0.305 

0.236 

3.0 

1.046 

0.908 

0.770 

10 

5.0 

1.469 

1.312 

1.151 

10.0 

2.888 

2.660 

2.471 

20.0 

6.5S5 

6.225 

5.572 

30.0 

11.014 

10.421 

9.175 

0 

0.802 

0.660 

0.423 

3.0 

1.460 

1.286 

0.946 

20 

5.0 

1.879 

1 .709 

1.319 

10.0 

3.271 

3.063 

2.563 

20.0 

6.911 

6.574 

5.667 

30.0 

11.264 

10.731 

9.208 

0 

1.104 

0.999 

0.525 

3.0 

1.747 

1.606 

1.033 

30 

5.0 

2.154 

2.009 

1.396 

10.0 

3.513 

3.342 

2.613 

20.0 

7.072 

6.S04 

5.659 

30.0 

11.430 

10.889 

9.141 

Table  13.  PITCHING  MOMENT  COEFFICIENTS 


( uearees » 

r— - - ■ - 

-30 


-20 


-10 


0 


10 


20 


30 


1 

MACH  2.60 

0 

7.313 

6.609 

3.437 

3.0 

7.308 

6.4S0 

3.601 

5.0 

7.016 

6.148 

3.451 

10.0 

6.632 

5.432 

3.254 

20.0 

6.353 

4.212 

2.740 

30.0 

2.980 

0.503 

1.222 

0 

5.307 

4.366 

2.797 

3.0 

5.185 

4.161 

2.828 

5.0 

4.7S2 

3.S02 

2.616 

10.0 

4.227 

3.152 

2.272 

20.0 

4.013 

2.317 

1.50S 

30.0 

0.745 

-1.259 

-0.188 

0 

2.550 

2.017 

m^QB 

3.0 

2.469 

1.944 

5.0 

2.146 

1.662 

HE;  tml 

10.0 

1.735 

1.160 

HUE  £K* 

20.0 

1.592 

0.414 

-0.034 

30.0 

-1.8S9 

-3.362 

-1.736 

0 

0.000 

0.000 

0.000 

3.0 

-0.043 

-0.034 

-0.054 

5.0 

-0.355 

-0.303 

BMW;' 

10.0 

-0.780 

-0.8IS 

20.0 

-1.167 

-1.S05 

-1.553 

30.0 

-5.033 

-6.229 

-3.093 

0 

-2.550 

-2.017 

-1.560 

3.0 

-2.655 

-2.115 

-1.586 

5.0 

-3.023 

-2.441 

-1.824 

10.0 

-3.646 

-3.150 

20.0 

-4.113 

-4.706 

-2.764 

30.0 

-7.570 

-8.942 

-3.960 

0 

-5.307 

-4.366 

3.0 

-5.398 

-4.61S 

ik Ok 

5.0 

-5.735 

-5.067 

10.0 

-6.1S2 

-5.816 

20.0 

-6.269 

-7.017 

-3.3SS 

30.0 

-9.224 

-10.995 

-4.178 

0 

-7.313 

-6.609 

-3.473 

3.0 

-7.297 

-6.732 

-3.326 

5.0 

-7.556 

-7.054 

-3.442 

10.0 

-7.784 

-7.664 

-3.479 

20.0 

-7.338 

-8.535 

-3.339 

30.0 

-10.323 

-12.038 

-3.732 

APPENDIX  B.  INTERPOLATED  WIND  TUNNEL  DATA 


Table  14.  AXIAL  FORCE  COEFFICIENTS  (POWER  ON) 


ALPHA 

(degrees) 

DESIGN 
POINT  1 

DESIGN 
POINT  2 

DESIGN- 
POINT  3 

El 

wmRmm 

0.280 

■■mi 

0.2S1 

0.284 

10.0 

0.296 

0.299 

0.336 

20.0 

0.335 

0.341 

0.367 

30.0 

0.373 

0.418 

0.349 

Table  15.  AXIAL  FORCE  COEFFICIENTS  (POWER  OFF) 


ALPHA 

(degrees) 

DESIGN- 
POINT  1 

DESIGN- 
POINT  2 

DESIGN 
POINT  3 

0 

0.326 

wmxsEmm 

3.0 

0.328 

5.0 

HU 

0.331 

10.0 

0.347 

■ MEM 

0.392 

20.0 

0.383 

mm 

0.421 

30.0 

0.421 

■91 

0.404 
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Table  16.  NORMAL  FORCE  COEFFICIENTS 


DELTA 

ALPHA 

DESIGN' 

DESIGN 

DESIGN 

(degrees) 

(degrees) 

POINT  1 

POINT  2 

POINT  3 

0 

iJtan 

-0.620 

-1.004 

3.0 

-0.090 

-0.395 

-30 

5.0 

0.2S5 

0.010 

10.0 

1.548 

1.362 

20.0 

4.813 

4.767 

4.883 

30.0 

8.712 

8.511 

9.013 

0 

-0.549 

Him 

-0.667 

3.0 

0.026 

-0.046 

-20 

5.0 

0.420 

Ilifil 

0.364 

10.0 

1.724 

1.707 

20.0 

5.052 

4.973 

5.173 

30.0 

8.953 

8.735 

9.284 

0 

-0.273 

-0.250 

-0.309 

3.0 

0.297 

0.299 

0.293 

-10 

5.0 

0.6SS 

0.6S5 

0.692 

10.0 

1.985 

1.969 

2.011 

20.0 

5.315 

5.217 

5.465 

30.0 

9.233 

8.987 

9.606 

0 

0.000 

0.000 

0.000 

3.0 

0.568 

0.550 

0.597 

0 

5.0 

0.95S 

0.936 

0.993  ! 

10.0 

2.259 

2.223 

2.316 

20.0 

5.603 

5.470 

5.S06 

30.0 

9.556 

9.235 

10.035 

0 

0.273 

0.309 

3.0 

0.S44 

0.915 

10 

5.0 

1.237 

■W 

1.320 

10.0 

2.547 

2.466 

2.671 

20.0 

5.920 

5.703 

6.243 

30.0 

9.839 

9.424 

10.451 

0 

0.549 

0.470 

0.66S 

3.0 

1.127 

1.014 

1.295 

20 

5.0 

1.527 

1.397 

1.718 

10.0 

2.S30 

2.663 

3.073 

20.0 

6.151 

5.848 

6.591 

30.0 

10.020 

9.513 

10.75S 

0 

0.77S 

0.620 

1.004 

3.0 

1.339 

1.148 

1.613 

30 

5.0 

1.723 

1.519 

2.016 

10.0 

3.002 

2.759 

3.351 

20.0 

6.26S 

5.88S 

6.817 

30.0 

10.073 

9.491 

10.916 
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Table  17.  PITCHING  MOMENT  COEFFICIENTS 


DELTA 

ALPHA 

DESIGN 

DESIGN 

DESIGN 

1 degrees) 

(degrees) 

POINT  1 

POINT  2 

POINT  3 

0 

5.145 

4.100 

6.644 

3.0 

5. 1 36 

4.177 

6.521 

-30 

5.0 

4.8S9 

3.990 

6.191 

10.0 

4.416 

3.690 

5.492 

20.0 

3.525 

3.036 

4.319 

30.0 

0.839 

1.078 

0.627 

|  | 

3.634 

3.111 

4.413 

3.539 

3.095 

4.212 

-20 

■IB 

3.24S 

2.853 

3.851 

10.0 

2.741 

2.448 

3.206 

20.0 

1.939 

1.670 

2.402 

30.0 

-0.759 

-0.402 

-1.159 

0 

1.SU4 

1.651 

2.044 

3.0 

1.748 

1.6ns 

1.9~0 

-10 

5.0 

1.4S0 

1.349 

1 .6S6 

10.0 

1.010 

0.903 

1.189 

20.0 

0.205 

0.056 

0.473 

30.0 

-2.603 

-2.061 

-3.2SS 

0 

0.000 

0.000 

3.0 

-0.050 

-0.034 

0 

5.0 

-0.309 

-0.313 

-0.306 

10.0 

-0.786 

-0.763 

-0.816 

20.0 

-1.6S7 

-1.603 

-1.773 

30.0 

-4.765 

-3.720 

-6.169 

0 

-1.804 

-1.651 

-2.044 

3.0 

-1.868 

-1.692 

-2.142 

10 

5.0 

-2.153 

-1.947 

-2.470 

10.0 

-2.669 

-2.377 

-3.175 

20.0 

-3.800 

-3.152 

-4.676 

30.0 

-6.617 

-4.956 

-8.873 

0 

-3.634 

-4.413 

3.0 

-3.746 

-4.657 

20 

5.0 

-4.072 

-5.100 

10.0 

-4.571 

-5.834 

20.0 

-5.323 

-4.050 

-6.980 

30.0 

-7.814 

-5.541 

-10.906 

-5.145 

-4.100 

-6.644 

1 

-5.142 

-4.007 

-6.760 

30 

-5.368 

-4.164 

-7.079  j 

10.0 

-5.711 

-4.316 

-7.670 

20.0 

-6.110 

-4.378 

-8.475 

30.0 

-8.162 

-5.393 

-11.952 
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APPENDIX  C.  LIFT  AND  DRAG  COEFFICIENT  DATA 


Table  18.  LIFT  COEFFICIENTS  (POWER  ON) 


DELTA 

(degrees) 

ALPHA 

(degrees) 

-30 

0 

3.0 

5.0 

10.0 

20.0 

30.0 

-20 

0 

3.0 

5.0 

1O.0 

20.0 

30.0 

-10 

0 

3.0 

5.0 

10.0 

20.0 

30.0 

0 

0 

3.0 

5.0 

10.0 

20.0 

30.0 

10 

0 

3.0 

5.0 

10.0 

20.0 

30.0 

20 

0 

3.0 

5.0 

1 0.0 

20.0 

30.0 

30 

0 

3.0 

5.0 

10.0 

20.0 

30.0 

DESIGN 
POINT  1 


DESIGN 
POINT  2 


-0.549 

0.011 

0.394 

1.646 

4.633 

7.567 


-0.273 

0.282 

0.661 

1.903 

4.880 

7.810 


0.000 

6.553 

0.930 

2.173 

5.151 

S.0S9 


0.273 

0.S2S 

1.208 

2.457 

5.448 

8.334 


0.549 

1.111 

1.497 

2.736 

5.665 

8.491 


0.778 

1.323 

1.692 

2.905 

5.775 

8.537 


2.665 

5.416 

8.010 


Table  19. 

DELTA 

(degrees) 


-30 


-20 


-10 


0 


10 


20 


30 


LIFT  COEFFICIENTS  (POWER  OFF) 


ALPHA 

DESIGN 

DESIGN 

DESIGN 

(degrees) 

POINT  1 

POINT  2 

POINT  3 

0 

-0.77S 

-0.620 

-1.004 

3.0 

-0.233 

-0.107 

-0.414 

5.0 

0.303 

0.255 

-0.023 

10.0 

1.38S 

1.464 

1.273 

20.0 

4.392 

4.347 

4.445 

30.0 

7.334 

7.139 

7.603 

0 

-0.549 

-0.470 

-0.668 

3.0 

0.009 

0.056 

-0.066 

5.0 

0.3S9 

0.425 

0.330 

10.0 

1.638 

1.650 

1.613 

20. 0 

4.616 

4.541 

4.717 

30.0 

7.543 

7.333 

7.838 

0 

-0.273 

-0.250 

-0.309 

3.0 

0.279 

0.281 

0.273 

5.0 

0.656 

0.654 

0.656 

10.0 

1.895 

1.S79 

1.912 

20. 0 

4.S63 

4.770 

4.991 

30.0 

7.7S6 

7.551 

S.117 

0 

0.000 

0.000 

0.000 

3.0 

0.550 

0.532 

0.577 

5.0 

0.925 

0.904 

0.956 

10.0 

2.164 

2.129 

2.213 

20.0 

5.134 

5.00S 

5.312 

3o.O 

8.065 

7.766 

8. 489 

0 

0.273 

0.250 

0.309 

3.0 

0.8265 

0.780 

0.S94 

5.0 

1.203 

1.150 

1.2S2 

10.0 

2.44S 

2.368 

2.562 

20.0 

5.432 

5.227 

5.723 

30.0 

8.310 

7.929 

8.S49 

0 

0.549 

0.470 

0.668 

3.0 

1.108 

0.995 

1.274 

5.0 

1.492 

1.363 

1.67S 

10.0 

2.727 

2.562 

2.958 

20.0 

5.649 

5.363 

6.050 

30.0 

8.467 

8.0075 

9.115 

0 

0.778 

0.620 

1.004 

3.0 

1.327 

1.129 

1.591 

5.0 

1.688 

1.484 

1.975 

10.0 

2.896 

2.657 

3.232 

20.0 

5.759 

5.401 

6.262 

30.0 

8.513 

7.987 

9.252 

Table  20.  DRAG  COEFFICIENTS  (POWER  ON) 


DELTA 

(degrees) 

ALPHA 

(degrees) 

0 

-30 

3.0 

5.0 

10.0 

20.0 

30. 0 

-20 

0 

3.0 

5.0 

10.0 

20.0 

30.0 

-10 

0 

3.0 

5.0 

10.0 

20.0 

30.0 

0 

0 

3.0 

5.0 

10.0 

20.0 

30.0 

0 

10 

3.0 

5.0 

10.0 

20.0 

30.0 

20 

0 

3.0 

5.0 

10.0 

20.0 

30.0 

30 

0 

3.0 

5.0 

10.0 

20.0 

30.0 

DESIGN 
POINT  1 


.27S 

.268 

.311 

.556 

.196 

.679 


DESIGN 
POINT  2 


0.2  SO 


0.591 

2.043 

4.800 


0.27S 

0.295 

0.342 

0.636 

2.133 

4.940 


0.27S 

0.309 

0.365 

0.684 

2.231 

5.101 


0.27S 

0.324 

0.390 

0.734 

2.340 

5.243 


0.27S 

0.339 

0.415 

0.783 

2.419 

5.333 


0.27S 

0.350 

0.432 

0.813 

2.459 

5.360 


0.2S0 

0.284 

0.323 

0.596 

2.021 

4.730 


0.280 

0.296 

0.340 

0.636 

2.105 

4.856 


0.2S0 

0.309 

0.364 

0.680 

2.191 

4.980 


0.2S0 

0.322 

0.386 

0.723 

2.271 

5.074 


0.280 

0.334 

0.405 

0.757 

2.321 

5.119 


0.280 

0.341 

0.415 

0.774 

2.334 

5.108 


DESIGN 
POINT  3 


0.319 

0.299 

0.323 

0.567 

2.015 

4.S09 


0.319 

0.317 

0.353 

2.114 

4.944 


0.319 

0.335 

0.3S2 

0.680 

2.214 

5.105 


0.319 

0.351 

0.408 

0.733 

2.331 

5.320 


0.319 

0.367 

0.437 

0.795 

2.4SO 

5.528 


0.319 

0.387 

0.472 

0.865 

2.599 

5.681 


0.319 

0.404 

0.497 

0.913 

2.676 

5.760 
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Table  21.  DRAG  COEFFICIENTS  (POWER  OFF) 


DELTA 

ALPHA 

DESIGN 

DESIGN 

DESIGN 

(degrees) 

(degrees) 

POINT  1 

POINT  2 

POINT  3 

0 

0.310 

■MM 

0.374 

3.0 

0.317 

0.355 

o 

ro 

• 

5.0 

0.360 

Ha 

0.378 

0.597 

0.610 

0.623 

20.0 

2.006 

1.994 

2.066 

30.0 

4.721 

4.657 

4.856 

0 

0.310 

0.326 

0.374 

3.0 

0.330 

0.331 

0.373 

-20 

5.0 

0.367 

0.369 

0.409 

10.0 

0.641 

0.642 

0.682 

20.0 

2. OSS 

2.065 

2.165 

30.0 

4.841 

4.769 

4.992 

0 

0.310 

0.326 

0.374 

3.0 

0.344 

0.343 

0.391 

-10 

5.0 

0.391 

0.389 

0.43S 

10.0 

0.686 

0.683 

0.735 

20.0 

2.17S 

2.148 

2.265 

30.0 

4.9S1 

4.S95 

5.153 

0 

0.310 

0.326 

0.374 

3.0 

0.358 

0.356 

0.407 

0 

5.0 

0.414 

0.411 

0.464 

10.0 

0.734 

0.727 

0.788 

20.0 

2.276 

2.235 

2.3S1 

30.0 

5.143 
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APPENDIX  D.  WIND  TUNNEL  GRAPHS 


Figure  19.  Lift  Coefficient  versus  Angle  of  Attack  at  Design  Point  One 
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Figure  20.  Lift  Coefficient  versus  Fin  Deflection  Angie  at  Design  Point  One 
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Figure  23.  Lift  Coefficient  versus  Angie  of  Attack  at  Design  Point  Two 
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Figure  27.  Lift  Coefficient  versus  Angle  of  Attack  at  Design  Point  Three 
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28.  Lift  Coefficient  versus  Fin  Deflection  Angle  at  Design  Point  Three 
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Figure  29. 


Pitching  Moment  Coefficient  versus  Angle  of  Attack  at  Design  Point 
Three 
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Figure  30.  Pitching  Moment  Coefficient  versus  Fin  Deflection  Angle  at  Design 
Point  Three 


APPENDIX  E.  SIMULATION  COMPUTER  PROGRAM 


?V*Vc***tV***V?*,;V‘iV*V?Vf***,sV*VfVfVr**Vr***Vf*VfVf*V)r***<Wf**,sV,;V*y?***Vr7!r* 
tV V? Vr •* *, V Vf * V y? iV ‘V Vc V? * V ‘V V? V? Vr -rV Vr * V i:  “V Vr -A- -,V Vr -,V ‘V '.V V? "V V?  Vr  Vr  y?  V? y?  V Vr V? * V? Vr y ,  ******** 


THESIS  RESEARCH 


***  DSL  SIMULATIONS  FOR  TOO  *** 

***  AUTOPILOT  CONFIGURATIONS  *** 

***  AT  THREE  DESIGN  POINTS  *** 

***  *** 

***  KENNETH  E.  COCKERHAM  *** 

***  CAPTAIN,  UNITED  STATES  ARMY  *** 

***  NAVAL  POSTGRADUATE  SCHOOL  *** 

»'iWnV'«V^V«V<V tVvV VciVm' VfvV~<V «V 


FIXED  INPUT 
RUN  1 

*'-****v?*yrVf ****** *DES IGN  POINT 

************ ***m t s g i t r  PARAMETERS******************* 

PARAM  I  =  720.0,  M  =  450.0,  DAMP  =  -22.4,  FLA  =  799276.0,  FLD 
PARAM  MPA  =  -168736.0,  MPD  =  -191368.0 


=  88808.0 


********FEEDBACK  GAINS  FOR  CONFIGURATION  ONE********* 

PARAM  KA1  =  -0.094,  KF1  =  27.6,  KBR1  =  0.752,  KST  =  26.96 

********FEEDBACK  GAINS  FOR  CONFIGURATION  TWO********* 

PARAM  KA2  =  1.0,  KF2  =  20.1,  KBR2  =  -0.51? 

"V 

******FEED  FORWARD  GAINS  FOR  CONFIGURATION  TWO******* 

PARAM  K1  =  -0.094,  K2  =  0.004 

*vyf  y?y*  *vyr*w»  y*y*  *w*ywvy?y»  y»  #v  #w*v*  «w**w?y»  »**  yr  <v  y-y-yj*  y?y*  ypyr  y»  v?«v  y?  y^yr/Vy*  yty* 


* 

* 


DERIVATIVE 

* 

* 

************ AUTO P I LOT  CONFIGURATION  ONE*************** 
* 

* 

R1  =  KA1  *  STEP(0.  0) 

XI  =  R1  -  X14  -  X12 

X2  =  INTGRL(0.  0 ,X1) 

X3  =  X2  -  XU  -  X13 
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X4 

= 

INTGRL(0.  0  ,X3) 

X5 

= 

X4 

DELTA 1 

= 

LIMIT(-0.  35,0.  35, X5) 

X6 

= 

DELTA 1  *  MPD/I 

ALFDD1 

= 

X6  +  X7  +  X8 

ALFD1 

= 

INTGRL(0.  OjALFDDl) 

ALF1 

= 

INTGRL( 0.0, ALFD1 ) 

X7 

= 

ALFD1  *  DAMP 

X8 

= 

ALF1  *  (MPA+MPD)/I 

X9 

= 

DELTA 1  *  FLD/M 

X10 

= 

ALF1  *  (FLA+FLD)/M 

ACCEL1 

= 

X9  +  X10 

XI 1 

KBR1  *  ALFD1 

X12 

= 

KST  *  Xll 

X13 

= 

KF1  *  X4 

X14 

= 

KA1  *  ACCEL 1 

************AUTOPI LOT  COST  I GURATI ON  TW 0*************** 


tic 


R2 

= 

KA2  *  STEP( 0.  0) 

Y1 

= 

R2  -  Y16 

Y2 

= 

K1  *  Y1 

Y3 

= 

INTGRL( 0.  0,Y2) 

Y4 

=: 

K2  *  Y1 

Y5 

= 

Y4  +  Y3 

Y6 

= 

Y5  -  Y14  -  Y15 

Y7 

= 

INTGRL(0.  0 ,Y6) 

Y8 

= 

Y7 

DELTA2 

= 

LIMIT(-0.  35,0.  35.Y8) 

Y9 

= 

DELTA2  *  MPD/I 

ALFDD2 

= 

Y9  +  Y10  +  Yll 

ALFD2 

= 

INTGRLCO.  0.ALFDD2) 

ALF2 

= 

INTGRL( 0.  0.ALFD2) 

Y10 

= 

ALF2  *  (MPA+MPD)/I 

Yll 

= 

ALFD2  *  DAMP 

Y12 

= 

ALF2  *  (FLA+FLD)/M 

Y13 

= 

DELTA2  *  FLD/M 

ACCEL2 

= 

Y12  +  Y13 

Y14 

= 

ALFD2  *  KBR2 

Y15 

= 

Y7  *  KF2 

Y16 

= 

KA2  *  ACCEL2 

* 

* 

* 

CONTROL  FINTIM  =2.0,  DELT  =  0.0001 


SAVE  0.001,  ACCEL1.ACCEL2 
END 
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* 

* 

-'-RUN  2 

****************X)£  s  I GN  PO I  NT  TWO********************* 

********* **»v***f j  XSSI LE  PARAMETERS******************** 

PARAM  I  =  693,  M  =  370,  DAMP  =  -11.9,  FLA  =  238350.0,  FLD  =  19080.0 
PARAM  MPA  =  -43545.0,  MPD  =  -54431.0 
********FEEDBACK  GAINS  FOR  CONFIGURATION  ONE********* 

PARAM  KA1  =  -0.845,  KF1  =  38.1,  KBR1  =  -4.44,  KST  =  8.98 
********feedback  GAINS  FOR  CONFIGURATION  TWO********* 

PARAM  KA2  =  1.0,  KF2  =  41.05,  KBR2  =  -7.62 
******FEED  FORWARD  GAINS  FOR  CONFIGURATION  TWO******* 

PARAM  K1  =  -0.845,  K2  =  -0.057 
***************************************************** 

END 

* 

* 

* 


*RUN  3 

***************DE SIGN  POINT  THRE E**********lV********* 
***************M  I S  S I  LE  PARAMETERS******************** 

PARAM  I  =  687,  M  =  365,  DAMP  =  -9.1,  FLA  =  99695.0,  FLD  =  10657.0 
PARAM  MPA  =  -26986.0,  MPD  =  -29679.0 

******** feedback  gains  for  configuration  one********* 

PARAM  KA1  =  -3.06,  KF1  =  40.9,  KBR1  =  -12.37,  KST  =  10.18 
Vc ***** Vc*fe E D B AC K  GAINS  FOR  CONFIGURATION  TWO********* 

PARAM  KA2  =  1.0,  KF2  =  51.17,  KBR2  =  -12.83 
******FEED  FORWARD  GAINS  FOR  CONFIGURATION  TWO******* 

PARAM  K1  =  -3.06,  K2  =  -0.352 
***************************************************** 


END 


GRAPH  (G1  ,DE=TEK618)  TIME,  ACCELl(RU=l,L0=-0.  2,SC=0.  2) 
LABEL  (Gl)  STEP  RESPONSE,  CONFIGURATION  1,  DESIGN  POINT  1 

jV 

GRAPH  (G2  ,DE=TEK618)  TIME,  ACCEL1(RU=2) 

LABEL  (G2)  STEP  RESPONSE,  CONFIGURATION  1,  DESIGN  POINT  2 


GRAPH  (G3  ,DE=TEK618)  TIME,  ACCEL1(RU=3) 

LABEL  (G3)  STEP  RESPONSE,  CONFIGURATION  1,  DESIGN  POINT  3 


GRAPH  (G4  ,DE*TEK618)  TIME,  ACCEL1(RU=1 , 2 , 3) 

LABEL  (G4)  STEP  RESPONSE,  CONFIGURATION  1,  DESIGN  POINTS  1, 
* 

GRAPH  (G5  ,DE=TEK618)  TIME,  ACCEL2(RU=1) 

LABEL  (G5)  STEP  RESPONSE,  CONFIGURATION  2,  DESIGN  POINT  1 

* 

GRAPH  (G6  ,DE=TEK618)  TIME,  ACCEL2(RU=2) 

LABEL  (G6)  STEP  RESPONSE,  CONFIGURATION  2,  DESIGN  POINT  2 
* 


GRAPH  (G7  ,DE=TEK618)  TIME,  ACCEL2(RU=3) 

LABEL  (G7)  STEP  RESPONSE,  CONFIGURATION  2,  DESIGN  POINT  3 


GRAPH  (G8  ,DE=TEK618)  TIME,  ACCEL2(RU=1,2,3) 


2,  &  3 


62 


LABEL  (G8)  STEP  RESPONSE,  CONFIGURATION  2, DESIGN  POINTS  1,  2,  &  3 
* 

GRAPH  (G9  ,DE=TEK618)  TIME,  ACCEL1(RU=1 ,LO=-0.  2,SC=0. 2) ,ACCEL2(RU=1) 
LABEL  (G9)  STEP  RESPONSE,  CONFIGURATIONS  1  &  2, DESIGN  POINT  1 

Vr 

GRAPH  ( G10  ,DE=TEK618)  TIME,  ACCEL1(RU=2) , ACCEL2(RU=2) 

LABEL  (G10)  STEP  RESPONSE,  CONFIGURATIONS  1  &  2,  DESIGN  POINT  2 

Vr 

GRAPH  (Gil  ,DE=TEK618)  TIME,  ACCEL1(RU=3) ,ACCEL2(RU=3) 

LABEL  (Gil)  STEP  RESPONSE,  CONFIGURATIONS  1  &  2,  DESIGN  POINT  3 

■jV 

END 

STOP 


P 
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